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We report a measurement of muon-neutrino disappearance in the T2K experiment. The 295-km 
muon-neutrino beam from Tokai to Kamioka is the first implementation of the off-axis technique 
in a long-baseline neutrino oscillation experiment. With data corresponding to 1.43xl0 20 protons 
on target, we observe 31 fully-contained single /i-like ring events in Super-Kamiokande, compared 
with an expectation of 104 ± 14 (syst) events without neutrino oscillations. The best-fit point for 
two-flavor — > v T oscillations is sin 2 (26 l 23) = 0.98 and |Am| 2 | = 2.65 x 10~ 3 eV 2 . The boundary 
of the 90% confidence region includes the points (sin 2 ^^), |Am| 2 |) = (1.0, 3.1x 10~ 3 eV 2 ), (0.84, 
2.65xlO~ 3 eV 2 ) and (1.0, 2.2x 10 _3 eV 2 ). 



PACS numbers: 14.60.Pq,13.15.+g,25.30.Pt,95.55.Vj 

We report a measurement of muon-neutrino disappear- 
ance in the T2K experiment. The muon-neutrino beam 
from Tokai to Kamioka is the first implementation of the 
off- axis technique [JJ in a long-baseline neutrino oscilla- 
tion experiment. The off-axis technique is used to pro- 
vide a narrow-band neutrino energy spectrum tuned to 
the value of L/E that maximizes the neutrino oscilla- 
tion effect due to Am^ , the mass splitting first observed 
in atmospheric neutrinos [2]. This narrow-band energy 
spectrum also provides a clean signature for subdomi- 
nant electron neutrino appearance, as we have recently 
reported [3]. Muon-neutrino disappearance depends on 
the survival probability, which, in the framework of two- 
flavor — > v T oscillations, is given by 

P surv = 1 - sin 2 (20 23 ) sin 2 (^f^) , (1) 

where E is the neutrino energy and L is the neutrino 
propagation length. We have neglected subleading os- 
cillation terms. In this paper we describe our observa- 
tion of disappearance, and we use the result to mea- 
sure |Am| 2 | and sin 2 (26*23). Previous measurements of 
these neutrino mixing parameters have been reported 
by K2K [lj and MINOS [5], which use on-axis neutrino 
beams, and Super-Kamiokande [BJ, which uses atmo- 
spheric neutrinos. 

Details of the T2K experimental setup are described 
elsewhere [7] . Here we briefly review the components rel- 
evant for the oscillation analysis. The J-PARC Main 
Ring (MR) accelerator [5] provides 30 GeV protons with 
a cycle of 0.3 Hz. Six bunches (Run 1) or eight bunches 
(Run 2) are extracted in a 5-/xs spill and are transported 
to the production target through an arc instrumented by 
superconducting magnets. The proton beam position, 
profile, timing and intensity are measured by 21 elec- 



trostatic beam position monitors (ESM), 19 segmented 
secondary emission monitors (SSEM), one optical tran- 
sition radiation monitor (OTR) and five current trans- 
formers. The secondary beamline, filled with helium 
at atmospheric pressure, is composed of the target, fo- 
cusing horns and decay tunnel. The graphite target is 
2.6 cm in diameter and 90 cm (1.9 \nt) long. Positively- 
charged particles exiting the target are focused into the 
96-m long decay tunnel by three magnetic horns pulsed at 
250 kA. Neutrinos are primarily produced in the decays 
of charged pions and kaons. A beam dump is located at 
the end of the tunnel and is followed by muon monitors 
measuring the beam direction of each spill. 

The neutrino beam is directed 2.5° off the axis between 
the target and the Super-Kamiokande (SK) far detector 
295 km away. This configuration produces a narrow-band 
beam with peak energy tuned to the first oscillation 
maximum E v = | Am§ 2 |L/(27r) ~ 0.6 GeV. 

The near detector complex (ND280) [7] is located 
280 m downstream from the target and hosts two de- 
tectors. The on-axis Interactive Neutrino GRID (IN- 
GRID) [9] records neutrino interactions with high statis- 
tics to monitor the beam intensity, direction and profile. 
It consists of 14 identical 7-ton modules composed of 
an iron-absorber/scintillator-tracker sandwich arranged 
in 10 m by 10 m crossed horizontal and vertical ar- 
rays centered on the beam. The off-axis detector re- 
constructs exclusive final states to study neutrino in- 
teractions and beam properties corresponding to those 
expected at the far detector. Embedded in the refur- 
bished UA1/NOMAD magnet (field strength 0.2 T), it 
consists of three large-volume time projection chambers 
(TPCs) [10] interleaved with two fine-grained tracking 
detectors (FGDs, each 1 ton). It also has a 7r°-optimized 
detector and a surrounding electromagnetic calorimeter. 
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The magnet yoke is instrumented as a side muon range 
detector. 

The SK water-Cherenkov far detector [TT] has a fidu- 
cial volume (FV) of 22.5 kt within its cylindrical inner 
detector (ID). Enclosing the ID is the 2 m-wide outer de- 
tector (OD). The front-end readout electronics [7j allow 
for a dead-time-free trigger. Spill timing information, 
synchronized by the Global Positioning System (GPS) 
with < 150 ns precision, is transferred from J-PARC to 
SK and triggers the recording of photomultiplier (PMT) 
hits within ±500 /is of the expected neutrino arrival time. 

The results presented in this Letter are based on the 
first two physics runs: Run 1 (Jan-Jun 2010) and Run 2 
(Nov 2010-Mar 2011). During this time period, the 
MR proton beam power was continually increased and 
reached 145 kW with 9 x 10 13 protons per pulse. The 
fraction of protons hitting the target was monitored by 
the ESM, SSEM and OTR and found to be greater than 
99% and stable in time. A total of 2,474,419 spills was 
retained for analysis after beam and far-detector qual- 
ity cuts, corresponding to 1.43 x 10 20 protons on target 
(POT). 

We present the study of events in the far detector 
with a single muon-like (/x-like) ring. The event selec- 
tion enhances charged-current quasi-elastic interac- 
tions (CCQE). For these events, neglecting the Fermi 
motion, the neutrino energy E v can be reconstructed as 



,xl0' 



E v = 



ml - (m» - Ebf 



+ 2(m„ - E b )En 



2(m n - E b - + p^ cos 6^ 



(2) 



where m p is the proton mass, m n the neutron mass, and 
Efj = 27 MeV the binding energy of a nucleon inside a 
16 nucleus. In Eq.J^ji^, and 6* M are respectively the 
measured muon energy, momentum and angle with re- 
spect to the incoming neutrino. The selection criteria for 
this analysis were fixed from Monte Carlo (MC) studies 
before the data were collected. The observed number of 
events and spectrum are compared with signal and back- 
ground expectations, which are based on neutrino flux 
and cross-section predictions and are corrected using an 
inclusive measurement in the off-axis near detector. 

Our predicted beam flux (Fig. [I]) is based on models 
tuned to experimental data. The most significant con- 
straint comes from NA61 measurements of pion produc- 
tion |12j in (p, 9) bins, where p is the pion momentum and 
the polar angle with respect to the proton beam; there 
are 5%-10% systematic and similar statistical uncertain- 
ties in most of the measured phase space. The production 
of pions in the target outside the NA61-measured phase 
space and all kaon production are modeled using FLUKA 
[l3l fT4] . The production rate of these pions is assigned 
systematic uncertainties of 50%, and kaon production un- 
certanties are estimated to be between 15% and 100% 
based on a comparison of FLUKA with data from Eichten 
et al. [H>]. The software package GEANT3 Hi . with 



o 
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1 2 3 4 5 

Neutrino Energy (GeV) 
FIG. 1. (Top) The predicted flux of z/ M as a function of neu- 
trino energy without oscillations at Super-Kamiokande and 
at the off-axis near detector; (Bottom) the flux of and 
at Super-Kamiokande. The shaded boxes indicate the total 
systematic uncertainty for each energy bin. 



GCALOR 17J for hadronic interactions, handles parti- 
cle propagation through the magnetic horns, target hall, 
decay volume and beam dump. Additional systematic 
errors in the neutrino fluxes are included for uncertain- 
ties in secondary nucleon production and total hadronic 
inelastic cross sections, uncertainties in the proton beam 
direction, spatial extent and angular divergence, the horn 
current, and the secondary beam line component align- 
ment uncertainties. The stability of the beam direction 
and neutrino rate per proton on target are monitored 
continuously with INGRID and are within the assigned 
systematic uncertainties [3]. 

Systematic uncertainties in the shape of the flux as 
a function of neutrino energy require knowledge of the 
correlations of the uncertainties in (p, 6) bins of hadron 
production. For the NA61 pion-production data [T^j, we 
assume full correlation between (p, 9) bins for each indi- 
vidual source of systematic uncertainty, except for par- 
ticle identification where there is a known momentum- 
dependent correlation. Where correlations of hadron- 
production uncertainties are unknown, we choose correla- 
tions in kinematic variables to maximize the uncertainty 
in the normalization of the predicted flux. 

Neutrino interactions are simulated using the NEUT 
event generator [18]. Uncertainties in cross sections of 
the exclusive neutrino processes are determined by com- 
parisons with recent measurements from the SciBooNE 
QU, MiniBooNE gDJUT], and K2K [H[23] experiments, 
comparisons with the GENIE [21] and NuWro [35] gen- 
erators and recent theoretical work |26) . 

An inclusive charged-current (CC) measurement in 
the off-axis near detector (ND) is used to constrain the 
expected event rate at the far detector. From a data sam- 
ple collected in Run 1 of 2.88 x 10 19 POT, neutrino inter- 
actions are selected in the FGDs with charged particles 
entering the downstream TPC. The most energetic neg- 
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atively charged particle in the TPC is required to have 
ionization energy loss compatible with that of a muon. 
The analysis selects 1529 data events with 38% CC 
efficiency and 90% purity. The agreement between the 
reconstructed neutrino energy in data and MC is shown 
in Fig. [2] The ratio of measured CC interactions to 
MC is 



R 



v^CC 
ND 



N 



Data,u ^CC 
ND 



N 



MCv^CC 
ND 



1.036 ± 0.028(stat.) 



±g;gf£(det.syst.) ± 0.038(phys.syst.), (3) 



where N ND 

T MC,v u CC 



Data,u^CC 



N 



ND 



is the number of i/n CC events, and 
is the MC prediction normalized by POT. 
The detector systematic errors in Eq.[3]are mainly due to 
uncertainties in tracking and particle identification effi- 
ciencies. The physics uncertainties result from cross sec- 
tion uncertainties but exclude normalization uncertain- 
ties that cancel in a far/near ratio. 



ED V, CCQE 
gg3 v„ CC non QE 

E3NC 
v„CC 
_J Outside FGD 




1.5 2 2.5 3 3.5 4 4.5 5 
Reconstructed V energy (GeV) 



FIG. 2. Neutrino energy reconstructed for the CCQE hy- 
pothesis for Vfj, CC candidates interacting in the FGD target. 
The data are shown using points with error bars (statistical 
only) and the MC predictions are in shaded histograms. 

At the far detector we select a CCQE enriched sam- 
ple. The SK event reconstruction [57] uses PMT hits in 
time with a neutrino spill. We select a fully-contained 
fiducial volume (FCFV) sample by requiring no activity 
in the OD, no pre-activity in the 100 [is before the event 
trigger time, at least 30 MeV electron-equivalent energy 
deposited in the ID, and a reconstructed event vertex in 
the fiducial region. The OD veto rejects events induced 
by neutrino interactions outside of the ID, and events 
where energy escapes from the ID. The visible energy re- 
quirement rejects events from radioactive decays in the 
detector. The fiducial vertex requirement rejects parti- 
cles entering from outside the ID. Further conditions are 
required to enrich the sample in CCQE events: a sin- 
gle Cherenkov ring identified as a muon, with momentum 
p p > 200 MeV/c, and no more than one delayed elec- 
tron. The muon momentum requirement rejects charged 
pions and misidentified electrons from the decay of un- 



TABLE I. Event reduction at the far detector. After each 
selection criterion is applied, the number of observed (Data) 
and MC expected events of i/ M CCQE, CC non-QE, intrin- 
sic u e , and neutral current (NC) are given. The columns de- 
noted by Ufj, include P M . All MC CC samples assume — > v T 
oscillations with sin 2 (20 2 3)=l.O and jAm 2 2 |=2.4 x 10~ 3 eV 2 . 

Data i^CCQE v^CC non-QE i/ e CC NC 



FV interaction 


n/a 


24.0 


43.7 


3.1 


71.0 


FCFV 


88 


19.0 


33.8 


3.0 


18.3 


single ring 


41 


17.9 


13.1 


1.9 


5.7 


/i-like 


33 


17.6 


12.4 


<0.1 


1.9 


p M > 200 MeV/c 


33 


17.5 


12.4 


<0.1 


1.9 


or 1 delayed e 


31 


17.3 


9.2 


<0.1 


1.8 



seen muons and pions, and the delayed-electron veto re- 
jects events with muons accompanied by unseen pions 
and muons. The number of events in data and MC after 
each selection criterion is shown in Table|T] The efficiency 
and purity of CCQE events are estimated to be 72% 
and 61% respectively. 

We calculate the expected number of signal events in 
the far detector (N^) by correcting the far-detector MC 



prediction with R^p from Eq. |ij 



Ns X K P (Er 



%c 2.^1 s ' 

E t 



(E t )N^(E r ,E t ). (4) 



In Eq. [ZJ Ng I I f(E ri E t ) is the expected number of events 
for the no-disappearance hypothesis for T2K Runs 1 and 
2 in bins of reconstructed (E r ) and true (Et) energies. 
Psurv(Et) is the two-flavor f M -survival probability, and is 
applied to and CC interactions but not to neutral- 
current interactions. 

The sources of systematic uncertainty in Ng X ^ are 
listed in Table HH Uncertainties in the near-detector and 
far-detector selection efficiencies are energy-independent 
except for the ring-counting efficiency. Uncertainty in 
the near-detector event rate is applied to N^^ a ' 1 '' 1 m 
Eq. [3j The flux normalization uncertainty is reduced be- 
cause of the near-detector constraint. The uncertainty in 
the flux shape is propagated using the covariance matrix 
when calculating Ng X ^ . The near-detector constraint also 
leads to partial cancellation in the uncertainty in cross 
section modeling, but the cancellation is not complete 
due to the different fluxes, different acceptances and dif- 
ferent nuclei in the near and far detectors. The total 
uncertainty in Ng X ^ is ^^'^ without oscillations and 

-14 8% witn oscillations with sm 2 (26> 23 ) = 1.0 and |Am| 2 | 
= 2.4° x 10" 3 eV 2 . 

We find the best-fit values of the oscillation param- 
eters using a binned likelihood-ratio method, in which 
sin 2 (2#23) and |Ato| 2 | are varied in the input to the cal- 
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TABLE II. Systematic uncertainties on the predicted number 
of SK selected events without oscillations and for oscillations 
with sin 2 (26> 23 ) = 1.0 and \Aml 2 \ = 2.4 x 10~ 3 eV 2 . 



Source 



(%, no osc) 



(%, with osc) 



SK CCQE efficiency 


±3.4 


±3.4 


SK CC non-QE efficiency 


±3.3 


±6.5 


orv 1NO efficiency 


1 n 


± I .A 


ND280 efficiency 


±5.5 -5.3 


±5.5 -5.3 


ND280 event rate 


±2.6 


±2.6 


Flux normalization (SK/ND280) 


±7.3 


±4.8 


CCQE cross section 


±4.1 


±2.5 


CCI71-/CCQE cross section 


±2.2 -1.9 


±0.4 -0.5 


Other CC/CCQE cross section 


±5.3 -4.7 


±4.1 -3.6 


NC/CCQE cross section 


±0.8 


±0.9 


Final-state interactions 


±3.2 


±5.9 


Total 


±13.3 -13.0 


±15.0 -14.8 



culation of 



until 



Nf£ a In 



\jdata 
AT exp 



(5) 



is minimized. The sum in Eq. [5] is over 50 MeV bins of 
reconstructed energy of selected events in the far detector 
from 0-10 GeV. 

Using the near-detector measurement and setting 
P surv = 1.0 in Eq. [IJ we expect a total of 103.6 t\tl 
(syst) single /u-like ring events in the far detector with- 
out disappearance, but we observe 31 events. If — > v T 
oscillations are assumed, the best-fit point determined 
using Eq. [| is sin 2 (26»2 3 ) = 0.98 and |Amf 2 | = 2.65 x 
10~ 3 eV 2 . We estimate the systematic uncertainty in 
the best-fit value of sin 2 (26*23) to be ±4.7% and that in 
I Ato§ 2 | to be ±4.5%. The reconstructed energy spectrum 
of the 31 data events is shown in Fig. [3] along with the 
expected far-detector spectra without disappearance and 
with best-fit oscillations. 

We construct confidence regions^] in the oscillation pa- 
rameters using the method of Feldman and Cousins |28j . 
Statistical variations are taken into account by Poisson 
fluctuations of toy MC datasets, and systematic uncer- 
tainties are incorporated using the method of Cousins 
and Highland [Ml ISO]- The 90% confidence region for 
sin 2 (2023) and |Am§ 2 | is shown in Fig. [d] for combined 
statistical and systematic uncertainties. 



In the T2K narrow-band beam, for a low-statistics data set, 
there is a possible degeneracy between the first oscillation maxi- 
mum and other oscillation maxima in L/E. Therefore we decided 
in advance to report confidence regions both with and without 
an explicit bound at |Am| 2 | < 5 X 10 — 3 eV 2 . For this data set, 
the bounded and unbounded confidence regions are identical. 



No-disappearance hypothesis 
Best-fit oscillation hypothesis 
T2K Run 1+2 data 




E 

z 



1 2 3 4 5 6 
Reconstructed energy (GeV) 
FIG. 3. Reconstructed energy spectrum of the 31 data events 
compared with the expected spectra in the far detector with- 
out disappearance and with best-fit z/ M — > v T oscillations. A 
variable binning scheme is used here for the purpose of il- 
lustration only; the actual analysis used equal-sized 50 MeV 
bins. 



,xl0" 3 



3.5 



> 



2.5 



1.5 



1 ' 

— • — T2K Run 1+2 




—• *.— T2K Run 1+2 Alternate Analysis 




MINOS 2011 




Super-K Zenith 




Super-K L/E 






. - V ~~ 















0.8 



0.85 



0.9 



0.95 



shT2e„ 



FIG. 4. The 90% confidence regions for sin 2 (26*23) and 
|Am| 2 |; results from the two analyses reported here are com- 
pared with those from MINOS [5] and Super-Kamiokande 

[snip. 



We also carried out an alternate analysis with a max- 
imum likelihood method. The likelihood is defined as: 



L 



Lnorm(sin 2 (20 23 ), Am 2 2 , f) 
L shape( sin2 ( 26, 23)> Am 32> f )£ syst (f), (6) 



where the first term is the Poisson probability for the ob- 
served number of events, and the second term is the un- 
binned likelihood for the reconstructed neutrino energy 
spectrum. The vector f represents parameters related to 
systematic uncertainties that have been allowed to vary 
in the fit to maximize the likelihood, and the last term 
in Eq. [6] is a multidimensional Gaussian probability for 
the systematic error parameters. The result is consistent 
with the analysis described earlier. The best-fit point for 
this alternate analysis is sin 2 (26 i 2 3) = 0.99 and |Ato 2 2 | 
= 2.63 x 10" 3 eV 2 . The 90% confidence region for the 
neutrino oscillation parameters is shown in Fig. [4j 

In conclusion, we have reported the first observation 



7 



of v ^ disappearance using detectors positioned off-axis 
in the beam of a long-baseline neutrino experiment. 
The values of the oscillation parameters sin 2 (2^23) and 
|Ar7i3 2 | obtained are consistent with those reported by 
MINOS [5] and Super-Kamiokande 
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